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Abstract. High redshift radio galaxies (HzRGs) provide unique targets for the study 
of the formation and evolution of massive galaxies and galaxy clusters at very high 
redshifts. We discuss how efficient HzRG samples are selected, the evidence for strong 
morphological evolution at near-infrared wavelengths, and for jet-induced star forma- 
tion in the z = 3.800 HzRG 4C41.17. 

1 Introduction 

Radio sources are surprisingly effective beacons for identifying galaxies at ex- 
tremely high redshifts. Optical/near-IR campaigns during the past few years by 
several groups have resulted in the discovery of more than 120 radio galaxies at 
z > 2, including 17 with z > 3, and 3 with z > 4. At low redshifts powerful radio 
galaxies are uniquely identified with massive ellipticals. If this is true also at 
high redshift, as seems reasonable given the surprisingly good Hubble K — z re- 
lation for radio galaxies at < z < 4.4, then we should be able to systematically 
study the evolution of massive elliptical galaxies over large lookback times using 
samples selected by their radio emission. While recently developed techniques of 
finding very distant star-forming galaxies (e.g. Steidel et al. 1996) are yielding 
substantial galaxy populations at z ~ 3, radio galaxy samples remain the best 
means of finding the most massive galaxies at z ~ 3, and even higher. 

Hierarchical galaxy formation scenarios suggest that these massive galaxies 
are assembled from smaller structures at relatively late cosmic epochs (e.g. Baron 
and White 1987). Observations of HzRGs may thus provide a unique opportunity 
to study the beginning of this process. Furthermore, when they have been found, 
HzRGs may also be used for 'color-dropout' searches of galaxy clusters around 
them, and begin cluster evolution studies at very high redshift. 

2 How To Find High Redshift Radio Galaxies 
2.1 Ultra Steep Spectrum Sources 

It has been known for many years that radio sources with very steep spectra 
are mostly associated with faint, distant galaxies (e.g. Tielens et al. 1979; Blu- 
menthal & Miley 1979). An example of a steep radio spectrum, for the powerful 
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nearby radio galaxy Cygnus A, is shown in figure 1. One can also see that the 
radio spectrum steepens at higher frequency. If all powerful radio galaxies have 
such curving spectra, due to synchrotron and/or inverse Compton losses, then 
with increasing redshift we would observe the steeper parts of their spectra. 
Analysis of the 3CR radio sample indeed showed such an expected a — z cor- 
relation (van Breugel and McCarthy 1990; Fig 2 shows an updated version). A 
similar result has recently been reported for a sample of ultra-luminous radio 
galaxies at high rcdshifts (Carilli et al. 1998). 
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Fig. 1. The radio spectral energy distribution of Cygnus A, using single dish mea- 
surements from the literature {e.g. Kellermann, Pauliny-Toth and Williams 1969, and 
others). Note the spectral steepening with frequency. 



Thus this ultra-steep spectrum (USS) 'red radio colors' selection technique 
may be used as a very effective tool for identifying very high redshift radio 
galaxies. The first comprehensive USS search for HzRGs was begun by Chambers 
and collaborators and resulted in the discovery of the z = 3.800 radio galaxy 
4C41.17 (Chambers, Miley and van Breugel 1990). For many years this was the 
most distant galaxy known. This relatively bright object, with its extended Ly- 
a halo, radio-aligned UV continuum, and evidence for star formation (see below) 
in many ways stimulated the continued search for more USS HzRGs. 

With the advent of several new, deep radio surveys it is now possible to define 
much better USS samples which reach 10 - 100 times fainter flux density lev- 
els than previous samples, and which allow much more accurate spectral index 
determinations. We have used several such surveys to define USS samples for 
study in the northern and southern hemispheres. Our primary northern hemi- 
sphere sample uses the WENSS 325 MHz survey (Rcngelink et al. 1997) together 
with the NVSS 1.4 GHz survey (Condon et al. 1998), to define the spectral in- 
dices, and the FIRST 1.4 GHz survey (Becker et al. 1995), to obtain radio maps. 
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Fig. 2. Spectral steepening with redshift for all 3CR radio galaxies. 



In the southern hemisphere we use the low frequency Texas 365 MHz (Douglas 
et al. 1996; Dec > -35°) and MRC 408 MHz (Large et al. 1981; Dec < -35°) 
surveys, in combination with the NVSS and PMN 5 GHz (Griffith & Wright 
1993; Dec < -35°) surveys. 

The sources in our southern hemisphere sample all have spectral indices 
a < —1.2. We have obtained radio images of all our southern USS targets with 
the Australia Telescope and VLA so that we can make accurate identifications. 
A deep near-infrared identification program of this sample was begun with the 
CTIO 4m telescope. Spectroscopic observations of some objects from this sample 
with the ESO 3.6 m have shown that USS sample is indeed extremely efficient 
at finding HzRGs, and have already resulted in the discovery of the most distant 
radio galaxy in the Southern Hemisphere known to date (TN J1338— 1942 at 
z = 4.13; De Breuck et al. , these proceedings). A summary of the current status 
of HzRG identification programs is given in Table 1 . 



Sample 


Definition 


Known Redshift 


Unknown 






z < 2 


2 < z < 3 


z > 3 




3CR 
MRC/Uy 
4C USS 
New USS 


Sits > 10 Jy 
S 4 os > 0.95 Jy 
alVs 4 < -1.0 
alf 5 ° < -1.3 


99.5 % 
93.4 % 
53 % 
35 % 


0.5 % 
5.9 % 

35 % 
35 % 


% 
0.7 % 
12 % 
30 % 


1% 
41 % 
50 % 
45% (faint) 



Table 1. The HzRG content for four samples of radio sources as a function of z: 3CR 
(Spinrad, private communication), MRC/Uy (McCarthy et al. 1996), 4C (Chambers 
et al. 1996), and our new USS sample. 
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Fig. 3. Hystogram of southern USS sources at DEC < -35°from the MRC and PMN 
surveys. In total 76 out of 2915 sources have atfgMflz < —1-20. 



2.2 The K — z Diagram 

Around the time that 4C41.17 was discovered another systematic search for 
HzRGs also proved to be successful. This method uses the magnitude - redshift 
relationship at infrared wavelengths, the K — z 'Hubble' diagram, for powerful 
radio galaxies galaxies to find promising targets. This resulted in the discovery of 
the z = 3.395 radio galaxy B2 0902+34 (Lilly 1988). Indeed, the K - z diagram 
has since proved to be a very reliable, if little understood, tool. Despite strong 
morphological evolution seen at near-infrared wavelengths, discussed below, the 
most powerful radio sources continue to follow the K — z relationship even to z <~ 
4.4 (Fig 4). It is our experience that, by combining the USS and K — z techniques 
and using near-infrared identifications of USS sources which are unidentified 
on UKST/POSS2 plates, one is virtually guaranteed to be successful in HzRG 
hunts. Without near-IR identifications and photometry approximately 2/3 of 
the galaxies are at z > 2. With near-IR identifications and photometry one can 
select the redshift range one wishes to study and to date all galaxies observed 
at Keck were found to be at the redshift predicted from the K — z diagram, at 
least for z < 4.4. 

3 Morphological evolution of HzRGs 

Near-infrared images obtained with the W. M. Keck I Telescope of HzRGs 
with 1.9 < z < 4.4 show strong morphological evolution at rest-frame opti- 
cal (A ros t > 4000 A) wavelengths (van Breugel et al. 1998; Fig 5). At the highest 
rcdshifts, z > 3, the rest-frame visual morphologies exhibit structure on at 
least two different scales: relatively bright, compact components with typical 
sizes of ~10 kpc surrounded by large-scale (~ 50 — 100 kpc) diffuse emission. 
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Fig. 4. Hubble K-z diagram for the 3C and 6C surveys (Eales et al. 1997), and NIRC 
data for HzRGs from van Breugel et al. (1998). The new USS A"-band detections and 
very high redshift candidates are shown by horizontal bars at the predicted redshift. 
The horizontal arrow indicates the range of predicted redshift for our faintest if-band 
object. The open star represents TN J1338-1942 (De Breuck et al. , this volume). All 
magnitudes are corrected to a 64 kpc metric diameter, assuming Ho = 50, qo = 0. The 
inset shows a blowup of the NIRC data in the 1.5 < z < 5.5 and 17.5 < K < 21 region. 

The brightest components are often aligned with the radio sources, and their 
individual luminosities are Mb ~ —20 to —22. For comparison, present-epoch 
L* galaxies and, perhaps more appropriately, ultra-luminous infrared starburst 
galaxies, have, on average, Mb ~ —21.0. The total, integrated rest-frame B- 
band luminosities are 3 — 5 magnitudes more luminous than present epoch L+ 
galaxies. 

At lower redshifts, z < 3, the rest-frame optical morphologies become smaller, 
more centrally concentrated, and less aligned with the radio structure. Galaxy 
surface brightness profiles for the z < 3 HzRGs are much steeper than those of 
at z > 3. We attempted to fit the z < 3 surface brightness profiles with a de Vau- 
couleurs r 1 / 4 law and with an exponential law, the forms commonly used to fit 
elliptical and spiral galaxy profiles, respectively. We demonstrate the fitting for 
our best resolved object at z < 3, 3C 257 at z = 2.474 (Fig 5). Within the limited 
dynamical range of the data, both functional forms fit the observed profiles — 
neither is preferred. Interestingly, despite this strong morphological evolution 
the K — z 'Hubble' diagram for the most luminous radio galaxies remains valid 
even at the highest redshifts, where a large fraction of the K-band continuum is 
due to a radio-aligned component. 
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Fig. 5. Selected near-IR images of HzRGs, presented in order of decreasing redshift, 
and the surface brightness profile of 3C257. 

4 Jet-induced star formation in 4C41.17: HST imaging 

Deep HST images of 4C41.17 at 'R'-band (rest-frame UV) and in Ly-a show that 
the 4C41.17 system consists of two components: 4C41.17-North with a bright 
string of UV knots and Ly-a emission along the radio axis, and 4C41.17-South 
with several much fainter UV knots, distributed in random fashion through- 
out a low surface brightness halo. The brightest radio knot in 4C41.17-North 
is associated with the brightest UV knot and arc-like Ly-a emission. One of 
the field objects in the HST images was also seen at near-IR (Graham et 
al. 1994; object # 16) and radio wavelengths (Carilli et al. 1994). This was used 
to align the HST and radio frames with an estimated relative accuracy of ~ 
0.1" . The central, radio-aligned UV and Ly-a emission is shown in Figure 6 
with the 0.21" resolution radio X-band image from Carilli et al. overlaid. Figure 
7 shows the HST continuum image of the entire 4C41.17 system smoothed to 
0.3" resolution to enhance low surface brightness features. 
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Fig. 6. HST WFPC2 images of 4C41.17-North with the X-band map of Carilli et 
al. 1994 superimposed. 



Downstream from the bright radio/UV/ Ly-a knot the radio source curves 
towards a faint, very steep spectrum NE lobe (see Carilli et al. 1994), while 
upstream from this knot, towards the central radio core, the UV continuum 
appears edge-brightened. These morphological features suggest a strong inter- 
action between the radio jet and dense ambient gas and, in fact, are as expected 
in jet-induced star formation models where sideways shocks induce star forma- 
tion in the dense medium of forming galaxies (e.g. De Young 1989; Bicknell et 
al. 1998). 

The star formation rates (SFR) of the various components as deduced from 
the rest-frame UV HST photometry range from 5-40 M©/yr for the knots 
in 4C41.17-North, to 5 - 10 M©/yr in 4C41.17-South. Here we have assumed 

L iwoi ~ 1040 1 ergs" 1 ^ 1 for a SFR = 1 M o/y r ( Conti et aL 1996 ) and 110 
dust reddening. The derived values are surely lower limits, given the detection 
of dust at sub — mm wavelengths in 4C41.17 by Dunlop et al. (1994). 

The entire 4C41.17 system is embedded in a common halo of diffuse, low 
surface brightness emission which extends over a very large area of 54/i^q 1 kpc x 
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76/igo 1 kpc (5" x 7" ). This includes a faint region, 4C41.17-South, with half a 
dozen compact knots distributed in random fashion. Spectroscopic observations 
have shown that 4C41.17-South is indeed at the same redshift as 4C41.17-North 
(Dey et al. 1999). The range of UV luminosities and SFR rates for the individual 
knots in 4C41.17-South is lower than in 4C41.17-North, and very similar to the 
'normal' (radio-quiet) Lyman-break galaxies discovered by Steidel et al. (1996). 
The random distribution and on average lower SFR in the 4C41.17-South knots 
suggests that star formation here is unaided by bowshocks from the radio jet. 
The total star formation rate, integrated over the entire 4C41.17 system and 
including the low surface brightness emission, is ~ 660 M Q /yr. Of this perhaps 
as much as 2/3 of the star formation may be occurring in the inter- knot regions. 
If the total star formation would continue at this rate for 2 x 10 8 — 2 x 10 9 yrs 
an entire massive elliptical galaxy of 10 11 M Q — 10 12 M Q might be assembled 
between z ~ 4 and z ~ 2.5, consistent with the morphological evolution for 
HzRGs seen in the near-IR Keck observations. 




Fig. 7. Smoothed version of the HST WFPC2/F702W image showing the entire 
4C41.17 system, including the clumpy companion system 4C41.17-South. 



5 Jet-Induced Star Formation In 4C41.17: Keck 
Spectroscopy 

Deep spectropolarimetric observations with the Keck II telescope by Dey et 
al. (1997) have provided strong evidence in support of the jet-induced star for- 
mation model for 4C41.17-North suggested above on the basis of the HST and 
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radio morphologies. These observations showed that the bright, radio-aligned 
rest-frame UV continuum is unpolarized (P;yy(2er) < 4%). This implies that 
scattered AGN light, which is generally the dominant contributor to the rest- 
frame UV emission in z ~ 1 radio galaxies, is unlikely to be a major component 
of the UV flux from 4C 41.17. Instead, the total light spectrum shows absorption 
lines and P-Cygni-like features that are similar to those detected in the spectra 
of the recently discovered population of star forming galaxies at slightly lower 
(z ~ 2 - 3) redshifts (Fig 8) . The detection of the S VA1502 stellar photospheric 
absorption line, the shape of the blue wing of the Si IV profile, the unpolarized 
continuum emission, the inability of other AGN-related processes to account for 
the UV continuum flux, and the overall similarity of the UV continuum spectra 
of 4C 41.17 and the nearby star forming region NGC 1741B strongly suggest 
that the light from 4C 41.17 is dominated by young, hot stars. The presence 
of radio-aligned features in many of the z > 3 HzRGs suggests, by analogy 
to 4C41.17, that jet-induced star formation may be a common phenomenon at 
these very high redshifts. 




Fig. 8. Keck spectrum (Dey et al. 1997) of the radio-aligned component 4C41.17-North, 
compared with a UV spectrum of the Wolf-Rayet starburst galaxy NGC 1741 (Conti 
et al. 1996). 
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